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SPECTROSCOPIC CONFIRMATION OF THREE z-DROPOUT GALAXIES AT z = 6.844 
DEMOGRAPHICS OF LYMAN-ALPHA EMISSION IN z - 7 GALAXIES* 
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ABSTRACT 

We present the results of our ultra-deep Keck/DEIMOS spectroscopy of z-dropout galaxies in the 
SDF and GOODS-N. For 3 out of 11 objects, we detect an emission line at ~ l/xm with a signal-to- noise 
ratio of ^ 10. The lines show asymmetric profiles with high weighted skewness values, consistent with 
being Lya, yielding redshifts of z = 7.213, 6.965, and 6.844. Specifically, we confirm the z — 7.213 
object in two independent DEIMOS runs with different spectroscopic configurations. The z = 6.965 
object is a known Lya emitter, IOK-1, for which our improved spectrum at a higher resolution yields 
a robust skewness measurement. The three z-dropouts have Lya fluxes of 3 x 10~^^ erg s~^ cm~^ 
and rest- frame equivalent widths EWq'*'" = 33 — 43 A. Based on the largest spectroscopic sample of 43 
z-dropouts that is the combination of our and previous data, we find that the fraction of Lya-emitting 
galaxies (EWg^" > 25A) is low at z ~ 7; 17 ± 10% and 24 ± 12% for bright (Muv ^ -21) and faint 
(Muv — —19.5) galaxies, respectively. The fractions of Lya-emitting galaxies drop from z ^ 6 to 7 
and the amplitude of the drop is larger for faint galaxies than for bright galaxies. These two pieces of 
evidence would indicate that the neutral hydrogen fraction of the IGM increases from z '^ 6 to 7, and 
that the rcionization proceeds from high- to low-density environments, as suggested by an inside-out 
reionization model. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: evolution — galaxies: 
high-redshift — 



1. INTRODUCTION 

Over the last two years, we have witnessed an explo- 
sion of activities aimed at searching for very high red- 
shift (z > 7) galaxies. This has been made possible 
by the installation of the Wide-Field Camer a 3 (WFC3) 
on-board the Hubble Space Telescope (e.g.. lOesch et al.l 
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Hawk-I instrument on the Very Large Telescope (VLT; 
iCastellano et al.|[2010ai rbl). and improvement in the red- 
end se nsitivity of the Su prime-Cam on the Subaru tele- 
scope (jOuchi et al.l 120091 ). These studies have identified 
over a hundred candidates at z > 6 .5 through the Ly- 
man break drop -out technique (e.g.. ISteidel et al1ll996l : 
IGiavaliscoll2002D , showing a decrease in the number den- 
sity of bright high-redshift galaxies w i th redshift (e.g., 
Ouchi et all 120091: iM cLurc ct al.' '2010^ 'Castcllano et alj 



2010aHBouwens et al . 2010a; Wilkins et al. 2010), bluer 



UV continua (e.g., [Bouwens et al.l l2010b[f^ . and rela- 
tively smaller stellar masses compared to l ower redshift 
galaxies selected by similar techiiiques (e.g.,lLabbe et al.l 
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Ono et al.l 120101) . Study of the intrinsic properties of 
high-redshift galaxies at epochs close to the dark ages is 
essential for understanding one of the most outstanding 
questions in modern astronomy - when did the Universe 
become reionized and what sources were responsible for 
it? This can be accomplished by studying the state of 
the inter-galactic medium (IGM) through estimates of 
the ionizing photon budget and Lyman a escape frac- 
tion. 

Several independent studies of the reionization pro- 
cess in recent years have yielded different results. By 
measuring the polarizati o n of the cosmic background 
radiation, IDunklev et al.l (|2009l) estimated the optical 

^^ See also lDunlop et aTl II2011I ). who question this result. 
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depth to reionization and concluded that, if it was a sud- 
den event, reionization occurred at z = 11.0 ±1.4 (see 
also. lKomatsu et al.ll2011l : lLarson et al.ll2011[ ). However 
investigating the spectra of SDSS quasars, iFan et aLl 
(|2006D studied the evolution of the Gunn-Peterson op- 
tical depth and demonstrated that the IGM reioniza- 
tion may have ended as la te as z ~ 6. The result from 
this study is questioned bv lGoto et al.l (|2011[ ). who have 
shown that for quasars at z > 6, a statistically large 
number is needed in Az = 0.1 bins to constrain cosmic 
variance and trace the evolution of the optical depth. 

Another useful tool for studying reionization is the 
Lyman a luminosity function of high-r edshift galaxies 
selected by narrow -band imaging (e.g., iHu et al.l 119981 : 
iRhoads et al.l I2000D : these galaxies are called Lyman 
Alpha Emitters (LAEs). Since neutral hydrogen in 
the IGM resonantly scatters Lya photons, the trans- 
mission of Lya is sensitive to the ionization state 
of the IGM. Therefore, we expect a decrease in the 
numbe r density of LAEs c l ose to t he reionization epoch 
e.g., iHai man fc Spaani 119991: iMalhotra fc Rhoadsl 
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2008; 
|2009l 120111) . 
decrease in the Lya 
luminosity function of LAEs, corresponding to an upper 
limit of the IGM neutral fraction xhi ^ 0.2 ± 0.2 at 
z = 6.6, which indicates that the major reionization 
process took place at hig her redshift. This resu lt has 
been further supported bv iKashikawa et ahl (|2011[ ). who 
studied the evolution of the L AE luminosity function 
in the Subaru Deep Field (SDF; IKashikawa et al. 2004) 
and found an increase in the neutral fraction o f the IGM 
from z = 5.7 to 6.5. iNakamura et al.l (|2011| ) have also 
found a large deficit in the number density of 2 = 6.5 
LAEs in the SSA22 field, and attributed it to significant 
ficld-to-ficld variance of the neutral fraction of the 
IGM. The effect o f cosmic variance was quantified by 
I Ouchi et al.l ()2010l ). who claimed a factor of 2 — 10 range 
in the number density of LAEs in an extensive survey 
covering an area of 1 deg^. 

Measuring the fraction of LAEs among Lyman-break 
Galaxies (LBGs; iGiavaliscol 120021) . the Lya fraction, 
provides complementary information to un derstand the 
reionization process (e.g.. lStark et al.ir2010[ ). Since LBGs 
are selected over a broader range of redshifts compared 
to observations with a narrow-band filter, their number 
density is less sensitive to cosmic variance. Searching for 
Lya emission fro m samples of LBGs with available spec- 
tra at 4 < z < 6, iStark et al.l ()2011f ) showed that lower 
luminosity LBGs have larger Lya fraction s, and that the 
fraction increases with redshift (see also. [ Vanzella et al.l 
I2OOI IStark et al.ll20Tol : iDouglas et al.ll2010( ). However, it 
is not yet clear if this trend continues z > 6. To explore 
this, we require spectroscopy of dropout candidates at 
z ~ 7. 

The real challenge in estimating the number density of 
z > 6 galaxies and the intensity of ionizing Lya photons 
is the spectroscopic confirmation of these candidates. 
They are extremely faint, and the only detectable fea- 
ture i s Lyg emission shifted to near-infrared wavelengths 
(e.g., Ilveet al.ll2006l : iLehnert et alll2010l : iFontana et all 



I2OIOI : IVanzella et"al]|2011f ). Despite significant efforts to 
spectroscopically confirm z > 6 candidat es, the number 
of con firmed sources is still very limited. IFontana et al.l 
(|2010f l reported the detection of one LBG with Lya emis- 
sion at z = 6.97 in ultra-deep spectroscopy of seven z- 
dropout candidates selected from VLT/Hawk-I imaging 
of the Great Observato ries Origins Deep Survey 's south- 
ern (GOODS-S) field (ICastellano et al.l I2010at) . They 
found a significant decline in the fraction of LBGs with 
Lya emission between z ^ 6 and 7, reversing the increas- 
ing trend with redsh ift found at z < 6. Furthermore, 
IVanzella et al.l (|2011l ) spectroscopically confirmed two z- 
dropout galaxies at z « 7.1 s elected from VLT/Hawk - 
I imaging of the BDF4 field (|Castellano et al.ll2010bf ). 
With the small number of z-dropout galaxies with avail- 
able spectroscopic redshifts, any measure of their Lya 
photon budget or number density will be seriously af- 
fected by statistical uncertainties and cosmic variance. 

In this study we present results from ultra-deep spec- 
troscopy with Keck/DEIMOS for a sample of 11 z- 
dropout galaxies. To maximize the spectroscopic success, 
we designed the photometric survey to identify dropout 
candidates at the bright end of the UV luminosity func- 
tion. The aim is to derive the fraction of LBGs with Lya 
emission and to study its evolution to z ~ 7. 

In the next section we present the photometric selec- 
tion of z-dropout candidates which are the targets for 
our spectroscopic observations here. The spectroscopic 
observations are described in Section |3l This is followed 
by redshift identification and the measurement of spec- 
troscopic properties in Section JH In Section |5| we discuss 
the implications for reionization. The conclusions are 
presented in Section [6l Throu ghout this paper, we use 
magnitudes in the AB system ()Oke fc Gu nn 198j) and 
assume a flat universe with (Sim, ^A, h) = (0.3, 0.7, 0.7). 

2. PHOTOMETRIC SELECTION OF 2-DROPOUT 
CANDIDATES 

In order to achieve successful spectroscopy of high red- 
shift LBGs, we need to identify a sample of bright candi- 
dates. However, the bright end of the luminosity function 
exponentially decreases. Therefore, we need to cover a 
large area to find a sufficient number of LBGs bright 
enough for spectroscopy. 

We carried out a wide-area photometric survey aimed 
at identifying z-dropout galaxies (i.e., galaxy candidates 
at z ~ 7) using Suprime-Cam on the Subaru telescope, 
outfitted with a custom-ma de filter (y-ha.nd) with ef- 
fective wavelength at l/im (jOuchi et al.l I2009D . This 
combination is ideal for identifying the bright popula- 
tion of LBGs at high redshifts. We covered an area of 
1568 arcmin^ to y ~ 26.0 mag (4(7 limit) for two fields: 
the SDF and the GO ODS northern field (GOODS-N; 
iGiavalisco et al.ll2004D . We identified 22 z-dropout can- 
didates with y = 25.4-26.1 mag, i.e., Muv < -21. This 
includes a galaxy which was alr eady identified to be at 
a spectroscopic redshift of 6.96 (jive et al.l 120061 ). These 
provide the targets for the spectroscopic observations in 
this paper. 

3. SPECTROSCOPIC OBSERVATIONS 

We used t he DEep Imaging Mult i- Object Spectrograph 
(DEIMOS; iFaber et all I2003D at the Nasmyth focus of 
the 10 m Keck II telescope to perform spectroscopic 
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TABLE 1 

Summary of observations with Keck/DEIMOS. 



Date (UT) 



Total Exposure 

fsccl 



Nj 



grating 
[lines mni^-'^] 



central wavelength 

[A] 



filter 



SDFZDIB SDF 2010 February 13 19350 5 830 

SDFZD3 SDF 2010 April 14 - 15 30000 2 830 

SDFZD4 SDF 2010 April 15 7200 4 830 

GNZDIB GOODS-N 2010 February 13, April 14-15 18000 2 830 

HDFllC GOODS-N 2011 April 1-2 14600 2 600 

HDFllD GOODS-N 2011 April 3 7200 2 830 



9000 
9000 
9000 
9000 
7500 
8100 



OG550 
OG550 
OG550 
OG550 
GG455 
OG550 



Numbers of observed 2-dropouts. Some objects were observed on multiple masks (See Tablc[2]|. 
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Fig. 1. — Spectra of SDF-63544 (left) and SDF-46975 (right). The top panels show the composite two-dimensional spectra, from which the 
one-dimensional spectra shown in the bottom panels are derived. A prominent emission line is seen at 9683 A for SDF-63544 {S/N ~ 13), 
and at 9536 A for SDF-46975 [S/N ~ 14). 



observations of th e z-dropout candidates discovered in 
lOuchi et all ()2009D . The data were taken on UT 2010 
February 13, April 14 - 15, and 2011 April 1-3. We 
observed 11 out of the 22 z-dropout candidates. We also 
observed the standard stars G191B2B and Wolf 1346 for 
flux calibration. The seeing was in the range 0.5" — 0.7". 
We used a total of six DEIMOS masks, as listed in Ta- 
ble [H For five masks, we used the OG550 filter and the 
830 lines mm~^ grating, which is blazed at 8640 A and 
was tilted to place a central wavelength of 9000 A on the 
detectors. This configuration provided a spectral cov- 
erage between 7000 A and 10400 A. For the remaining 
mask (HDFllC) we used the GG455 filter and the 600 
lines mm~^ grating, blazed at 7400A. This was tilted 
to place a central wavelength of 7500A on the detector. 
The spectral coverage was between 5200 A and 10200 A. 
The spatial pixel scale was 0.1185" pix~^, and the spec- 
tral dispersion was 0.47 A pix~^ and 0.65 A pix~^ for 
the 830 lines mm~^ and 600 lines mm~^ grating, respec- 
tively. The slit widths were 1". For objects filling the 
slit, the FWHM resolution of the 830 and 600 grating 
was ~ 3.3A and ~ 4.7A, respectiveljO- Details of the 

^* http : //www2 . keck . hawaii . edu/inst/deimos/specs . html 



observations, the filters, gratings and the total exposure 
time used arc listed in Table [TJ 

We used the GG455 filter for HDFllC, although this 
filter allows transmission of second-order light redward of 
~ 9100 A. This filter was used because some of the tar- 
gets on those masks were z ^ 3 ultra-luminous infrared 
galaxy candidates, which would have spectral signatures 
below 5000 A. Although this configuration allowed con- 
tamination from second-order light of an emission line 
whose wavelength is longer than ~ 4500A, the targeted 
z-dropouts in the HDFllC mask were also observed in 
other configurations without such contaminations in the 
red wavelength range. 

The reduction was performed using the spec2d IDL 
pipelin j^ developed by the DEEP2 Redshift Survey 
Team (jDavis et al.ll2003|) . We used a modified version 
of the spec2d (Capak et al. in prep.) for the data taken 
with HDFllD, since those data were dithered. Wave- 
length calibration was achieved by fitting to the arc lamp 
emission lines. The spectra were flux calibrated with the 
standard stars G191B2B and Wolf 1346. We applied no 

^^ The pipeline was developed at UC Berkeley with 
support from NSF grant AST-0071048. Downloaded at 

http : //astro . berkeley . edu/~cooper/deep/spec2d/ 
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Fig. 2. — Spectrum of the z = 7.213 z-dropout galaxy, GN- 
108036. The top panels show its two-dimensional spectra obtained 
with the GNZDIB, HDFllC, and HDFllD masks. The size along 
the spatial axis is 5.0" for each two-dimensional spectrum. The 
HDFllD spectrum is binned in 2 X 2 pixels. A line is visually 
identified at ~ 9980A in the spectra of GNZDIB and HDFllC, 
whose exposure times are 5 hours and ~ 4 hours, respectively, 
while the line is marginally seen in the spectrum of HDFllD, whose 
exposure time is 2 hours. In the bottom panel, we show the one- 
dimensional spectra. The gray solid lines are spectra obtained with 
individual masks. The composite spectrum is shown as the black 
solid line. All the one-dimensional spectra illustrate a line detection 
at around 9980 A, and the S/N ratio of the line in the composite 
spectrum is ~ 6. 

correction for slit-loss effects, since the same slit width 
was used for both the standard stars and science targets. 
We estimate the uncertainty due to the slit-loss correc- 
tion by comparing the fluxes of the two-dimensional spec- 
tra in the slit width with total fluxes and find this to be 
less than 10%. We measure the la sky noises from the 
pixel distributions around 9600 A, the wavelength corre- 
sponding to that of Lya at z « 6.9, the expected peak 
of the r edshift distribution for z-dropout selection (Fig- 
ure 6 of lOuchi et al.l[2009| ). In the measurements, we do 
not avoid the wavelength ranges significantly affected by 
strong OH lines. We find the la sky noise to be (1.4—5.7) 
xlO^^^ erg s~^ cm~^ (for details, see Section |4!2|) . We 
summarize the la flux limits in Tabled 

4. RESULTS 
4.1. Redshift Identification 

Figure [T] shows a prominent emission line in the spectra 
of SDF-63544 and SDF-46975. We fit Gaus sian profiles 
to these lines using the IDL MPFIT routine (jMarkwardtl 
l2009fP^ . and find the emission lines correspond to the 
observed central wavelengths 9683A (SDF-63544) and 
9536A (SDF-469750- One z-dropout galaxy, SDF- 

^® http://www.physics.wisc.edu/~craigm/idl/ldl.html 
^"^ The best-fit central wavelengths are slightly shorter (~ 3 A), 
if we fit an asymmetric Gaussian profile in which the standard 
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Fig. 3. — Two-dimensional spectrum of GN-108036 in the wave- 
length range of 7400 — 7800A. The size along the spatial axis is 
5.0". No emission line is detected in this wavelength range, which 
argue against the possibility that the detected line at ~ 9980A is 
Ha at 2 = 0.521. In this case, an [Olll]A5007 line would have been 
detected at approximately 7617A marked by an arrow (b). No line 
detection also suggests that the detected line at ~ 9980A is not 
either an [Olll]A5007 at z = 0.994 or H/3 at z = 1.054. In the 
former/latter case, the [Oll]AA3726, 3729 doublet would have been 
detected at the position marked by arrows (a)/(c). 

63544, was previously identified as a Lya emitter, lOK- 
1, whose redshift was confirmed spectroscopically by 
live et all (J2006I ). However, the signal-to-noise ratio 
{S/N) of the Lya line in their spectrum is only ~ 5. 
We deemed it worthwhile to detect the line at higher 
S/N and with higher resolution in order to verify the 
line detection and identify the line based on line profile 
analysis. We confirm that the central wavelength of the 
emission line is almost the same as that derived previ- 
ously (9682A). 

For one of the targets in GOODS-N, GN-108036, we 
detect an emission line at about 9980 A (Figure [2]). We 
confirm the line detection in three independent DEIMOS 
observations, performed in 2010 and 2011, with different 
configurations. In 2010, we obtained the spectrum us- 
ing mask GNZDIB, the 830 lines mm~^ grating, and the 
OG550 filter (first two-dimensional spectrum in the top 
panel in Figure [2). In the 2011 run, we used a different 
set up with the 600 lines mm~^ grating, the GG455 filter, 
and a different mask (HDFllC) to locate the spectrum 
on a different position on the DEIMOS CCD (second one 
in the top panel in Figure [2]). The HDFllD data were 
obtained over two nights of observing, and the line is 
detected in independent reductions of each night's data, 
as well as in the combined spectrum. We also used a 
third configuration for our 2011 run, with the 830 lines 
mm~^ grating, the OG550 filter, and an entirely new 
mask HDFllD (third one in the top panel in Figure [2]). 
Although the grating and the filter here are the same as 
those used in the 2010 observation, the mask is differ- 
ent. The 2011 observations with the 830 grating were 
also dithered, whereas the 2010 observations with the 
830 grating were undithered, and the pipeline processing 
was correspondingly different for the two data sets. The 
three observations of GN-108036 independently confirm 
that the line is detected at three different positions on 
the DEIMOS CCDs, with two different spectroscopic se- 
tups. This strongly argues against the observed feature 
being an artifact. We fit Gaussian functions to the line 
profile in Figure [2l using the MPFIT routine, and find a 
central wavelength of 9984A. 

We investigate the possibility that the observed lines 

deviation on the blue side is smaller than that on the red side, 
as a better approximation of the observed profiles but with an 
additional free parameter. 
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TABLE 2 
Summary of Follow-up Spectroscopy 



Object 


Mask 


Total Exposure 


flux limit (Icr) 


y(total) 


EWn^" limit {3a) 






[sec] 


[erg s-i cm-2] 


[mag] 


[A] 


GN-152505 


GNZDIB, HDFllC, HDFllD 


39800 


3.0 X 10"^** 


25.2 


10 


GN-108036t 


GNZDIB, HDFllC, HDFllD 


39800 


4.2 X 10-1** 


25.5 


7.1 


SDF-63544t 


SDFZDIB, SDFZD3 


49350 


2.1 X 10-1* 


25.1 


5.8 


SDF-83878 


SDFZDIB 


19350 


3.3 X 10-" 


25.1 


12 


SDF-46975t 


SDFZDIB, SDFZD3 


49350 


1.9 X 10-1* 


25.2 


7.1 


SDF-76507 


SDFZDIB 


19350 


2.4 X 10-1* 


25.2 


9.2 


SDF-123919 


SDFZD4 


7200 


5.3 X 10-1* 


25.4 


31 


SDF-75298 


SDFZDIB 


19350 


1.4 X 10-1* 


25.5 


6.5 


SDF-121418 


SDFZD4 


7200 


4.5 X 10-1* 


25.6 


31 


SDF-107344 


SDFZD4 


7200 


5.7 X 10-1* 


25.7 


53 


SDF-136726 


SDFZD4 


7200 


5.1 X 10-1* 


25.7 


45 



Note. — The flux limits of the objects with Lya detections (GN-108036, SDF-63544, and SDF-46975) were 
estimated using the sky-noise distribution in the vicinity of the Lya line, while those of objects lacking Lya 
detections were estimated using the sky-noise distribution between 9400 — 9800 A. The EWq^" limits of the 
Lya-detected objects were estimated by dividing their Lya flux limits by their UV continuum flux densities and 
(1 -I- Zgpec), while those of objects without Lya detection were estimated assuming their redshifts are equal to 6.9. 

t 2:-dropout with Lya detection. 



TABLE 3 

Properties of Spec:;troscopically Confirmed LY'a-EMiTTiNG z-dropouts 



Object 



Rcdshift 



^Lyc 



[erg s- 



mcont 

[mag] 



[mag 



EW, 



Lya 



5™ FWHM Vfwhm 

[A] [A] [A] [kms-1] 



GN-108036 7.213 
SDF-63544 6.965 
SDF-46975 6.844 



2.5 X lO-i'^ 1.5 X 10*3 
2.8 X lO-i'^ 1.6 X 10*3 
2.7 X 10-1" 1.5 X 10*3 



25.2 
25.4 
25.4 



-21.8 
-21.6 
-21.5 



33 
43 
43 



4.1 
12.6 
8.6 



+0.7 
0.7 

+0.4 
0.4 

+0.6 
0.7 



15 

6.7 
12 



442 
207 
374 



are something other than Lya (i.e., Ha, H/3, [On], or 
[Oiii]). The main argument in favor of the observed 
lines being Lya is their morphology and the elear asym- 
metry of the lines (see below). Furthermore, the lines 
are unlikely to be H/3 or [Oiil] a.t z ^ 0.9 — 1.0, since, 
in this case, we expect to detect additional lines. If 
the detected line were H/3, the [Oiii]A5007 line would 
fall at 9972A for SDF-63544, 9822A for SDF-46975, and 
10284A for GN-108036. However, none of the objects 
show the corresponding detections, implying the line ra- 
tio [Oiii]A5007/H/3 of < 0.2 - 0.5 (Scr upper limit). 
Galaxies with 12 + log(0/H) > 8.8 - 9.0 meet these 
low ratios (e.g., Nagao et al. 2006), but our objects 
are unlikely to be such metal rich, because if they were, 
they should be very massive (from the mass-metallicity 
relation) and thus their broadband magnitudes would 
be much brighter than observed. In the case that the 
detected line were [Oiii]A5007, the H/3 line would be 
seen at 940lA for SDF-63544, 9259A for SDF-46975, 
and 9674A for GN-108036. However, none of the ob- 
jects show the corresponding detections, implying the 
line ratio [Oili]A5007/H/3 of > 2.0 - 4.0 (3cr lower limit). 
These lower limits correspond to su bsolar oxygen abun - 
dances (~ 7.0 - 8.7; See Figure 17 of lNagao et all[2006h . 
which may be typical of low-mass galaxies. Thus, un- 
fortunately the non-detection of H/3 cannot strongly rule 
out the possibility of our objects being [Oiii] A5007 emit- 
ters. In addition, the lines are unlikely to be [On]. If 
the lines were [On] emitters at z '^ 1.6, the FWHIVI reso- 
lution of 3.3A would have distinguished the two compo- 



ne nts of the double t, separated by ~ 7A. In their Figure 
10, iHu et al] ()2004D showed the spectra of emission-line 
objects with the [On] doublet signature, which were ob- 
tained with Keek/DEEVIOS using the same configuration 
as ours (the 830 grating and the OG550 filter). The fact 
that the detected lines in our spectra are singlet strongly 
argues against the possibility that they are z ~ 1.6 [On] 
emitters. Moreover, the non-detection of the galaxies in 
the deep i-band image^3 strongly disfavors the possi- 
bility that these galaxies are actually at z ^ 1.6. For 
GN-108036, if the detected line were Ha at z = 0.521, 
then we might also expect to detect [Oiii] or H/3, which 
is actually not detected as shown in Figure [H For SDF- 
63544 and SDF-46975, we cannot cheek the detection of 
[Oiii]A5007 since their spectra do not cover the wave- 
length range blueward of 7500A. 

To quantify the asymmetry of the lines, we intro- 
duce the weight e d ske wness parameter, S^, following 
iKashikawa et al.l (|2006l ). This can be used to distin- 
guish Lya from other emission lines. S^, is defined as 
the product of the skewness (the third moment of fiux 
distribution) and the width of the line. Lya emission 
lines at high redshifts typically have large positive S^ 
values, while other possible lines. Ha, H/?, and [Oili], 
are nearly symmetric and have almost zero values of 
Sw When the [On] line is not resolved, its 5*^, value 
is expected to be small, sin ce the [Oii]A3726 line is 
weaker than [Oii]A3729 (e.g.. iRhoads et al.l[200l . The 

1* The 2(T limiting total magnitude of the SDF (GOODS-N) 
i-band image is 28.3 (27.6). 
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Sw values of the three z-dropouts are estimated to be 



O-i/ 



12.6 ± O.4A for SDF-63544, S„ 



8.6 



+0.6 



A for 



SDF-46975, and 5*^ = 4.1 ± 0.7A for GN-108036, which 
means that all the three lines have an asymmetric pro- 
file with a sharp decline on the blue side and a long tail 
on the red sid e, as is commonly seen in Lyg at high 
redshifts (e.g., iShimasaku et al.l 120061 : iKashikawa et al.l 
|2006[) . IKashikawa et al.l (I2006D empirically set S^ = 3 A 
as a critical value to distinguish Lya emission from 
other emission lines f or galaxies at z > 5.7 (see also, 
IShimasaku et al.ll2006( ) . Our objects have higher S^ than 
the criterion, which would suggest that the detected lines 
are Lya. 

Note that the asymmetric line profile in the spectra 
of GN-108036 might be caused by over-subtractions of 
OH lines blueward of the detected line. However, there 
is a window of OH lines near the detected line at the 
blue side. If the line had a symmetric profile, then we 
would expect to see faint emission in this window, sim- 
ilar to that seen on the red side of the line profile, but 
none is actually detected. The sky subtraction might 
affect the line profile to some extent, but the asymmet- 
ric line profile is likely to be real. In addition, there is 
no line detected other than the one at ~ 9980A. Fur- 
thermore, the galaxy is not detected in the deep optical 
broadband images from either Subaru or HST/ACS. Its 
colors {z — y> 1.6 {2a), j/ — rnFi4ow = 0.3) are consistent 
with those expected for a LEG at z = 7.2. All of these 
factors support the interpretation that the detected line 
is redshifted Lya. 

Thus, we conclude that all of the detected lines are 
redshifted Lya. The redshifts derived from the Lya 



emission centroids are 



^spcc 



6.965 for SDF-63544, 



z,pe^ = 6.844 for SDF-46975, and Zspec = 7.213 for GN- 
108036. These redshifts might be overestimated since 
the Lya emission lines of LBGs are t ypically shifted red- 
ward of their systemic redshifts (e.g., l Pettini et al1l2001l : 



Liloii: 



IShaplev et ani2003l: ISteidel et al. l l2010h . The velocity off- 
sets are typically less than ^ 1000 km s~^, which cor- 
responds to a difference between the systemic redshift 
(zsys) and that estimated from the Lya line (^Lya) of 

Az = Zsys - Zi^ya < -0.027. 

Nv A1240 is the only high-ionization metal line, indica- 
tive of AGN activity, which would fall within our spec- 
tral range. The line would fah at 9876A for SDF-63544, 
9727A for SDF-46975, and 10184A for GN-108036. None 
of the objects shows this emission line, placing a 3cr lower 
limit to the line ratio Lya/Nv of > 4 for SDF-63544, 
> 5 for SDF-46975, and > 2 for GN-108036. However, 
since a typical high-z AGN has a line ratio of Lya/Nv 
= 4-20 (e.g., lMcGarthvlll993f ). these lower limits are 
mostly not useful to exclude the possibility that the ob- 
jects are AGN hosts; to do so much deeper spectroscopy 
is needed. Throughout this paper, we assume that the 
light of our objects is not contaminated by AGNs. 

4.2. Line Flux Measurement 

We compute the Lya line fiux by summing flux den- 
sities around the line, neglecting the minor contribu- 
tion of UV continuum. This corresponds to the range 
9677 - 9699 A for SDF-63544, 9526 - 9551 A for SDF- 
46975, and 9977-9997 A for GN-108036. The estimated 
fluxes arc summarized in Table |3l 




Fig. 4. — Images of GN-108036 taken by the Subaru Suprime- 
Cam (y), and HST WFC3 (F140W). The size of each panel is 
5"x5". North is up and east is to the left. 
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Fig. 5.— Images of GN-108036 in the IRAC 3.6iim (top) and 
4.5/im (bottom). The panels from left to right show the original 
IRAC images, the best-fit model images constructed from PSF 
templates at the object positions measured from the HST WFC3 
image, and the residuals after subtracting the model from the data. 
The size of each panel is 20.4" x20. 4". North is up and east is to 
the left. 



SDF -63544 was previously discovered by live et al.] 
(|2006[ ) who report Lya fluxes of 2.0 x 10^^^ erg s~^ cm~^ 
and 2.7 x 10"^'' erg s~^ cm~^ from their spectrum and 
narrow-band image, respectively. Our estimated Lya 
flux of SDF-63544 is 2.8 x 10"^'^ erg s^^ cm^^^ consistent 
with those previous measurements. 

In order to estimate flux limits for the objects with- 
out Lya detections, for each object we sample the one- 
dimensional spectra in 25A bins, comparable to the 
width of Lya lines. We then estimate the la flux limit 
by fltting Gaussian functions to the flux distribution over 
the wavelength range 9400 - 9800A. The la flux limits 
for individual objects are listed in Table [2] 

4.3. Equivalent Width Measurement 

We do not directly detect continuum in the spectra 
of the z-dropout galaxies. We can only estimate their 
continuum flux densities from their y-band magnitudes 
m easured over a 1.8 arcscc diameter aperture (Table 3 
of lOuchi et al.|[2009() . allowing for aperture corrections, 
IGM absorption and Lya emission. 

To estimate the aperture correction for our z-dropouts, 
we select bright, non-saturated and isolated sources in 
the y-band images, and perform multi-aperture photom- 
etry to construct the curve of growth. We then determine 
the aperture correction by measuring the difference be- 
tween magnitude over the 1.8 arcsec diameter aperture 
and the total magnitude (the asymptote of the growth 
curve). The median aperture correction is 0.344 mag for 
GOODS-N and 0.349 mag for the SDF. 
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TABLE 4 
SED Fitting Results for z = 7.213 galaxy 
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log{Age) 

[yr] 
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Fig. 6.— Results of SED fitting for the z = 7.213 galaxy GN- 
f08036. The filled squares show the observed flux densities used 
for SED fitting, while the open squares show those not used. The 
vertical arrows sho w 5(t upper limits for HST/ACS Vgoei *775i -^850 
llUuchi et al.|[2009r) . Subaru/MOIRCS J, H, and K from left to 
right. The left panel is for 'stellar + nebular' {f^°^ = 0) models. 
The red solid curve indicate the best-fit SED that is the sum of 
a stellar SED (red dotted curve) and a nebular SED (red dashed 
curve). The crosses indicate synthesized flux densities in individual 



bandpasses. The right panel is for 'pure stellar' (/^' 



1) models. 



The flux in the j/-band is affected by Lya emission 
and IGM absorption shortward of rest-frame 1216A. We 
correct for the contribution due to Lya emission using 
the estimated Lya flux from section H?^ for objects with 
strong Lya emission, and the 3cr flux limit for the z- 
dropouts without Lya detection. The IGM absorption is 
taken into account o n the assumpti on of flat UV contin- 
uum in /^ and using iMadaul (|1995D . The 3<t EW limits 
for the 11 galaxies in the spectroscopic sample are listed 
in Table [21 In Table [31 we list our estimates for the con- 
tinuum magnitude nicont at the Lya wavelength, and the 
corresponding EWs that we derived for the Lya emission 
lines. 

4.4. Photometry and Stellar Population Properties of 
the z = 7.213 galaxy 

The Subaru y-band photometry for the z = 7.2 galaxy 
GN-108036 is difficult to interpret because the flux mea- 
surement is strongly affected both by the Lya emission 
line and the Lya forest IGM absorption, which should 
suppresses roughly 65% of the intrinsic continuum flux 
from the galaxy within the y bandpass. Fortunately, GN- 
108036 has been imaged with the HST WFC3-IR camera 
(Weiner et al. in preparation; HST program #11600). 
The image has a total exposure time of 1217 sec and 
was obtained with the F140W filter, whose central wave- 
length is approximately lAjim, sampling 1700 A in the 
rest-frame UV continuum. The Suprime-Cam y-band 
and WFC3 F140W images of GN-108036 are shown in 
Figure [HI Photometry in the HST image measures a to- 
tal magnitude of 25.17 ± 0.07 mag. GN-108036 was also 
observed in the Subaru MOIRCS JH K survey of the 
GOODS-N field (jKaiisawa et al.|[20Tll ). The galaxy is 



not significantly detect in any of these MOIRCS bands, 
with 5(7 photometric upper limits of 24.9 (J), 24.4 (77), 
and 24.6 (A'), marginally consistent with the WFC3 pho- 
tometry. 

The WFC3 F140W magnitude of GN-108036, m == 
25.17, translates to a rest-frame 1700A luminosity 
-^1700 = —21.81 at z = 7.213. There have been many 
recen t estimates of the UV luminosity function at z = 7 
(e.g. lOuchi et al. 2009; ' Oesch et al.|[20Tol: iMcLure et all 



l2010t iGrazian et al.. .20111) . These have consistently 



found values of the characteristic luminosity M* in the 
range —19.9 to —20.3. GN-108036 is therefore an impres- 
sively luminous galaxy, with Luv ~ 4 to 6L*; indeed, it 
is roughly twice a s bright as an L*j^f LBG at z = 2 to 
3 (Af{Jv « -21.0. [Reddv et al.l l200'8). The exponential 
cut-off of a Schechter luminosity function would imply 
that galaxies this luminous would be rare indeed. For 
example, the best-fit Schechter parameterization of the 
z = 7 luminosity function from IBouwens et al.l (I2011b[ ) 
yields a space density 2.7 x 10~^ Mpc~'^ for galaxies with 
M\]Y < —21.81. For a redshift interval Az = 1 at z = 7, 
we would expect to find only 0.1 galaxies this luminous 
within the 160 arcmin^ GOODS-N ACS/WFC3/IRAC 
field, or about 1 galaxy over the 1568 arcmin^ covered 
by our Suprime-Cam survey in the combined SDF and 
GOODS-N fields. In fact, the redshift selection function 
for z-dropouts in o ur z — y survey is sig nificantly narrower 
than Az := 1 (see lOuchi et al.ll2009f ). making the large 
UV luminosity of GN-108036 still more remarkable. It is 
impossible to judge based on a posteriori statistics from 
a single object, but either we were quite lucky to find 
and spectroscopically confirm such a bright galaxy, or 
perhaps the bright end of the UV luminosity function at 
z = 7 has been underestimated in studies to date. In fact, 
the observational constraints on the bright end of the lu- 
minosity function come fr om just a few ground-based sur- 
veys, including our own ([Ouchi e t al. 200§) and that of 
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iCastellano et al.l ()2010al fbl). Evidently more deep imag- 
ing and spectroscopic studies over wider fields are needed 
to provide better measurements. The y-band itself be- 
comes an unreliable indicator of luminosity for galaxies 
at z > 7 without exact measurements of the galaxy red- 
shift and Lya line flux, so deep near-infrared data at 
longer wavelengths (such as that from HST WFC3) will 
be essential for a robust determination of the luminosity 
function. 

GN-108036 also falls within the field covered by the 
extremely deep IRAC imaging from the GOODS Spitzer 
Legacy program ( PI: M. Dickinson). As was noted in 
lOuchi et al.l ()2009l ). the galaxy is faintly detected in the 
IRAC 3.6/im and 4.5^m images, which roughly sample 
the rest-frame B and V band light at z = 7.2, although 
it is partially blended with a foreground galaxy located 
1" to the east (Figure [5]). In order to extract reliable 
IRAC fluxes, we measure the positions of the two galax- 
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ies in the WFC3 F140W image, and then position unit- 
normahzed IRAC PSF images at these locations in the 
background-subtracted 3.6/im and 4.5/ini images. We 
then fit these PSF templates to the IRAC data, minimiz- 
ing ^ to derive the best-fitting fiuxes and their uncer- 
tainties. The middle panels in Figure [5] show the best-fit 
model and the right panels are the residual images. In 
this way, we obtain the total magnitudes of GN-108036 
of ma.e^m = 25.44 ± 0.13 and m^.^^,^ = 24.86 ± 0.13, 
respectively. 

We fit models to the photometry of GN-108036 to in- 
fer its stellar po pulation properti es. The procedure is the 
same as that of lOno et al.l ()2010[ ). except that redshift is 
fixed at z = 7.213. We use the stellar population synthe- 
sis model GALAXEV (iBruzual fc Charlotl[200l for the 
stellar component of the SEDs, and consider two extreme 
cases for nebular emission: /p°" = 0, where ionizing pho- 
tons are all converted into nebular emission (the 'stellar 
-I- nebular' case), and Z^™ = 1, where all ionizing photons 
escape from the galaxy (the 'pure stellar' case). Nebu- 
lar spectra (lines and continua) are c alculated following 
the procedure given in lSchaerer fc de Barros (2009). We 
adop t the Salpeter initial mass function (IMF: iSalpeteii 
Il955| ) , a constant rate of star formation, and stellar and 
gas metallicities Z = 0.2 Z(^. For dust attenu ation, we 
use the functional form of iCalzetti et al.l (j20Q0f ) with the 
ass umption that E{B — F)gas = E{B — F)*, as proposed 
by lErb et al.l (|2006f ). The model SEDs are fit to the 
observed flux densities in the WFC3 F140W and IRAC 
3.6/im and 4.5/um bands, and its MOIRCS JHK magni- 
tudes. We do not use the y-band photometry since it is 
strongly affected by IGM attenuation and Lya emission. 
The free parameters in the fitting are stellar mass, age, 
and dust extinction, with three degrees of freedom. 

The results of SED fitting are summarized in Table lU 
and the best-fitting SEDs are shown in Figure |6l The 
best-fit models have small stellar masses, 4 x YQ^Mq to 
3 X 1O^M0, very young ages, 4 to 32 Mvr, and modest 
color excesses E{B - V)^ = 0.02 to O.IO Deeper K- 
band photometry would be needed to reliably constrain 
the amplitude of the Balmer break, but the blue rest- 
frame UV-to-optical color tofi4ow — "1,3. g^m = —0.27 
suggests that the break is small and that the light from 
the galaxy is dominated by very young and relatively 
unreddened stars. 

If the fiux in the IRAC bands has a significant contri- 
bution from nebular emission, then the intrinsic stellar 
SED would be even bluer than the directly measured 
WFC3-to-IRAC color would suggest (see Figure HI left). 
There is evidence that this is the case. The very red 
IRAC 3.6/im to 4.5/im color is not easily reproduced by 
stellar emission alone, but can be explained by a strong 
contribution of nebular emission in the 4.5/im bandpass. 
At z = 7.2, the IRAC 3.6/im band is mostly free of 
strong emission lines, but the IRAC 4.5/im band includes 
[Oiii]A4959 and A5007 (H^ is largely excluded). Several 
recent studies have suggested that strong nebular emis- 
sion lines can significantly affect IRAC photometry for 



high redshift galaxi e s (e.g.. iSchaerer fc de Barros 
iVanzella et al.lboiOt iRaiter et al.ll2010t lOno e%~. 



2009 : 



20101: 



^^ The uncertainties should be larger if star-formation history 
and metallicity are varied, although we fix them due to the small 
number of the data points. 



iShim et"aLl 120111) , and galaxies with extremely strong 
[Qui] (EWq > lOOO A) have been ident ified at z < 1 (e.g., 
iKakazu et aIll2007D and z = 1 to 1.5 (jAtek et al.l [20ia). 
Here we see evidence for strong [Oiii] at z = 7.2; in the 
best-fit stellar -|-nebular model, the [Om] lines contribute 
- 60% of the flux in the IRAC 4.5/im band. 

Using the conversion from UV c ontinuum lurn i nosity 
density to star formation rate from iMadau et ahl ()1998[ ) 
for a Salpeter IMF, and assuming no extinction, the 
I700A luminosity of GN-108036 implies a star forma- 
tion rate (SFR) of 29M0 yr~^. However, the best-fitting 
stellar population models give larger SFR k, IOOM0 yr""'^ 
(Table S]). In part this is due to extinction in the models, 
and in part to the very young ages that are implied by 
the SED-fitting. The standard co nversion factors of U V 
luminosity per u nit SFR (e.g., from lMadau et al.]ll998l or 
lKennicut3ll998[ ) assume star formation timescales > 10* 
years. At younger ages, the UV continuum is still build- 
ing toward its steady-state value, and SFRs derived using 
the standard conversion factors will be underestimated. 
For comparison, the Lya line luminosity (1.5 x 10**^ erg 
s~^) yields another estimate of SFR(LyQ;) = I3.6M0 
yr~^, again adopting a standard Salpeter IMF conversion 
factor (jKennicuttll 19981 ) and assuming Case B recombina- 
tion for the Ha/Lya flux ratio. This SFR is significantly 
smaller than either estimate from the continuum SED, 
as expected since Lya is subject to strong attenuation 
from dust extinction and by the IGM. 

Because the unusual IRAC color of GN-108036 
strongly favors the results from the 'stellar-|-nebular' 
model fitting, we conclude that the stellar mass of 
this galaxy is most likely smaller than lO^Af0. This 
value is not atypical compared t o other published esti- 
mates for galaxies at z ^ 7 (e.g.. lFinkelstein et al.|[2010l : 
[Gonzalez et al. 201l|). However, when combined with the 
bright UV luminosity and the large derived star forma- 
tion rate, it implies an extremely high specific star forma- 
tion rate (SSFR) > 10""^ yr"!. This is at least 50 times 
higher than the mean values tha t have been estimate d 
for samples of LBGs at 4 < z < 7 (jGonzalez et al.| [2010;l. 
and suggests that GN-108036 is seen at a special mo- 
ment in its life cycle when it is undergoing a very strong 
starburst, presumably with a short duration. 

Although the IRAC color favors 'stellar -|-nebular' mod- 
eling, the very young age derived may be in conflict with 
a simple expectation that the duration of star forma- 
tion cannot be shorter than the dynamical time of the 
system, since it will take around the dynamical time to 
have the gas of the system collapse to form stars. GN- 
108036 has an FWHM of 0.30" in the F140W image 
and a Lya line width of 4.4 x 10^ km s~^. Dividing 
the former by the latter gives a rough estimate of the 
dynamical time of '^ 4 Myr, which is marginally con- 
sistent with the age from SED fitting. Although this 
possible problem is not directly related to the discussion 
and conclusions of this paper, similar conflicts between 
stella r age and dynamical time are found in other paper s 
(e.g.. ISchaerer fc de Barr"oi[200l [2OIOI: lOno et al.l[2010l) . 
highlighting the need for further research. 

Assuming a flat /^ UV spectral slope, we derive an 
equivalent width for Lya EWq^" = 33 A. Another galaxy 
from our sample, SDF-63544, was also observed with 
HST/WFC3 in the F125W, F140N, and F160W band- 
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Fig. 7. — Examples of mock spectra made by artificially adding a mock emission line to the two-dimensional spectrum of GN-152505. 
The mock emission line has Lyo equivalent width of 50A and continuum magnitude of 25.57. The left panel shows a spectrum whose mock 
line is found at ~ 9660A, while a mock line {~ 9380A) is not recovered in the right spectrum due to the residuals of severe OH lines around 
the artificial line. 



passes ()Cai et al.l 1201 IL Egami et al. in preparation). 
A fit to this photometry yi elds a UV continuum slope 
/a oc A~^-^^ (jCai et al.ll2011f ). and a continuum flux mag- 
nitude of 25.5±0.1 at 9683A, the wavelength of Lya, con- 
sistent with our estimate from the j/-band photometry in 
Section 1431 



4.5. Lya Fraction 

In order to infer the ionizing state of the IGM, wc study 
the evolution of the Lya fraction by measuring its value 
at z ~ 7 and comparing it with previous meas urements at 
4 < z < 6. For 4 < 2 < 6 dropout galaxies. iStark et"al] 
(|201lD divided their sample into two UV luminosity bins, 
-21.75 < Muv < -20.25 and -20.25 < Afuv < -18.75, 
and estimated their Lya fractions. Since the UV absolute 
magnitudes of our z-dropouts arc in the brighter range, 
in the following analysis for our sample we focus only on 
the Lya fraction of dropout galaxies at —21.75 < Muv < 
-20.25. 

We define the Lya fraction as the ratio of the number 
of z-dropout galaxies with strong Lya emission measured 
from their spectra to the total number of spectroscopi- 
cally observed z-dropouts. 



vEWa Y^ Lya / \^ 



(1) 



where a, ^^ is 1 if any i-th galaxy has a Lya EW larger 
than a critical EW (EWc) and otherwise, pi is the prob- 
ability that the Lya wavelength of the i-th z-dropout 
candidate is within the observable wavelength range - 
i.e., the range not contaminated by OH airglow lines, 
and is given hy pi ~ J C[{z)dz/ J C{z)dz, where C{z) is 
the probability distribution for z-dropouts as a function 
of redshift and Cl{z) is the effective redshift probability 
distribution function of z-dropout candidates when cor- 
rected for OH emission. Such a statistical analysis of the 
Pi parameter is needed since, at the wavelength range 
~ 9000 — 10100 A, where we expect to detect their Lya 
line, our z ^ 7 candidates may not always satisfy the 
EWc criterion (i.e., pi = 1), as the atmospheric OH air- 
glow lines will significantly affect the spectrum in that 
wavelength range. 

In order to compute C'{z), we perform Monte Carlo 
simulations by generating mock Lya emission lines with 
EWo = 50A or lOOA and y-band magnitudes of 25.24 
and 25.57 mag. These magnitudes correspond to the 
central values of two bins into which the y-band total 
magnitudes of our z-dropouts are divided in order of 
their brightness. We add the mock Lya lines to an ob- 
served two-dimensional spectru m, following the redshift 
distribution, C{z) (Figure 6 in lOuchi etahl |2009l Wc 
then inspect the mock two-dimensional spectrum search- 



ing for a Lya feature. Figure [7] shows examples of two- 
dimensional spectra apparently with and without Lya 
detected (left and right plots, respectively). We perform 
a number of trials and evaluate the recovery rate of mock 
Lya lines. The C'{z) is then estimated by computing 
the weighted mean o f Lya EWs, assum ing this to be 
a Gaussian function (|Ouchi et al.ll2008() . We find that 
over 90% of the mock Lya lines are successfully recov- 
ered in our extremely deep spectra (i.e., C^{z) ~ C{z)). 
Therefore, we choose not to correct for this effect in the 
following analysis. When computing the Lya emission 
fraction, Xj^ " above a certain equivalent width thresh- 
old EWc, we divide the number of galaxies with detected 
EW(Lya) >EWc by the number of galaxies for which the 
3cr upper limits to the EW detectability are smaller than 
EWc (i-G-, those galaxies for which Lya lines stronger 
than EWc could have been detected). 

From the sample of 11 z-dropout galaxies with spec- 
troscopic data studied here, we detect Lya emission lines 
for 3 sources, all with Lya EWs larger than 25A. How- 
ever, we do not include GN-108036 at z = 7.2, since it 
is brighter than —21.75 in the rest-frame UV. We find 
4 objects with Lya EW limits larger than this, giving a 
Lya fraction, Xl^^ = 33 ± 27% (2/6). If we set the 

EWc = 55A, the 2 objects with Lya detection have 
smaller EWq^" than the criterion. The Lya EW upper 
limits of 8 z-dropout candidates without Lya detection 
arc lower than the EWc. We thus obtain an upper limit 
of Lya fraction: Xl^^ < 10% (0/10). 

5. DISCUSSION: IMPLICATIONS FOR REIONIZATION 

Since neutral hydrogen in the IGM resonantly scat- 
ters Lya photons, the transmission of Lya photons is 
sensitive to the ionization state of the IGM. Thus, the 
fraction of Lya-emitting LBGs can be used as a di- 
agnostic for the ionization state of the IGM. In Sec- 
tion 14.51 we estimated the Lya fraction of bright z- 
dropout galaxies to be X^^^ = 33 ± 27% (2/6), and 

^Lya < 10^0 (0/10)- In order to study evolution of the 
Lya fraction to z '^ 7, we compare our estimates for 
the bright (-21.75 < il/uv < -20.25, or Muv - -21) 
z-dropout galaxies with those at z ^ (4, 5, 6), Xl^^ = 
(12 ±3%, 24 ±5%, 20±8%) andXfj^^ = (5.5±1.8%, 
7.2±2.8%, 7.5±5.0%) and their extrapolations to z ~ 7, 
Xl^^ = 33 ± 10% and Xl^^ = 9 ± 6%, which are de- 
rived on the assumption of a linear relationsh ip between 
Lya fraction and redshift (jStark et al.l 120 111) . For faint 
(-20.25 < Muv < -18.75, o r Mtt v ~ -19.5) galaxies 
at z - (4, 5, 6), IStark et all ()201lh obtained their Lya 
fractions of Xl^^ = (35 ± 5%, 48 ± 9%, 54 ± 11%) and 
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Properties of 



TABLE 5 

-DROPOUTS SpECTROSCOPICALLY OBSERVED IN THE PREVIOUS STUDIES 



Object 



Rcdshift^' 



[mag] 



Muv" 
[mag] 



[erg s^-^ 



Ewf;^° 



comments 



[A] 



Fontaria ct al. (2010) 



G2_1408 
G2.2370 
G2_4034 
G2_6173 
H.9136 

W.6 

0.5 



BDF-521 
BDF-3299 



6.972 



26.37 


-20.49 


3.4 X IQ- 


re— 


25.56 


-21.27 


< 2.5 X IQ- 


-IM 


26.35 


-20.50 


< 2.5 X IQ- 


-IM 


26.53 


-20.33 


< 2.5 X IQ- 


-IM 


25.90 


-20.94 


< 2.5 X 10" 


-1« 


26.93 


-20.38 


< 2.5 X 10- 


-1« 


27.52 


-19.67 


< 2.5 X lo- 


-18 



7.008 
7.109 



25.86 
26.15 



-20.63 
-20.56 



Vanzella et ^ (2011)' 



13 S/N{Lya) = 7. Castcllano ct al. (2010a') . YbPEN Hawk-I 

< 4.8 Castellano et al. (2010a '). 'Kopen Hawk-I 

< 9.7 Castellano et al. (2010a ). 'Kopen Hawk-I 

< 11 ^astellano_^t_^l^ (2010a), Y^PEN Hawk-I 

< 6.4 Hickev et al. (2010') . 'Kopen Hawk-I 

< 11 'Wilkins et al. (2010) . Yogs 'WFC3-ERS 

< 21 Oesch et al. (2010). Yipl •WFC3-HUDF 



1.62 X 10"'^ 
1.21 X 10"^^ 



64 
50 



S/N(Lya) 
S/N{Lya) 



18 
16 



ERS5847 

ERS7376 

ERS7412 

ERS8119 

ERS8290 

ERS8496 

ERS10270 

ERS10373 

A1703jD1 

A1703jD3 

A1703jD6 

A1703jD7 

A2261.1 

BoRG.4 

EGS.K1 

HUDF09.799 

HUDF09.1584 

HUDF09.1596 

MS0451-03.10 



Schenker et al. 



6.48f 
6.79t 
6.38t 
6.78t 
6.52t 
6.441 
7.02t 
6.44t 
6.75+ 
6.89+ 
7.045 
8.80+ 
7.81+ 
8.27+ 
8.27+ 
6.88+ 
7.17+ 
6.905 
7.50+ 



26.6 
27.0 
27.0 
27.1 
27.1 
27.3 
27.4 
27.4 
24.1 
25.5 
25.8 
26.8 
26.9 
25.8 
25.3 
27.7 
26.7 
26.8 
26.7 



-20.22 
-19.89 
-19.79 
-19.79 
-19.73 
-19.51 
-19.54 
-19.41 
-20.39 
-19.25 
-19.36 
-18.74 
-18.85 
-21.39 
-21.89 
-19.21 
-20.27 
-20.12 
-16.10 



9.1 X 10" 



2.8 X 10" 



65 



65 



30 



Pcntcricci ct al. (20lTy" 



with a close neighbor 
S/N(Lya) > 5 



Magnification factor /x = 9.0+ 
Magnification factor /x = 7.3+ 
Magnification factor /^ = 5.2+, S/N{Lya) > 5 
Magnification factor fi = 5.0+ 
Magnification factor fj, = 3.5 



Magnification factor fi = 50 



3.2 X 10"^*^ 16~ 



NTTDF-474 

NTTDF-1479 

NTTDF-1632 

NTTDF-1917 

NTTDF-2916 

NTTDF-6345 

NTTDF-6543 

BDF4-2687 

BDF4-2883 

BDF4-5583 

BDF4-5665 



6.623 



6.701 



26.50 
26.12 
26.44 
26.32 
26.64 
25.46 
25.75 
26.15 
26.15 
26.65 
26.64 



-20.35 
-20.77 
-20.45 

-20.25 
-21.41 
-21.14 
-20.74 
-20.74 
-20.24 
-20.25 



7.2 X 10" 



15 



S/NiLya) = 7 

likely an interloper 
S/N(Lya) = 11 



Note. — The upper limits of Lya flux and equivalent width are bcr. 
^ The redshift of objects without spec-z is set at 6.8. 
*= For the objects of lFontana et al.l | |20TO|). I'Vanzella et al.l i20Tii) . and IPentericci et al.l ll20Tll 



_ _j __ _ y-band magnitudes arc taken from their 

Table 1, respec tively. For the object s of lSchenker et al.l 1)20111 ). Ji7s magnitudes are taken from their Table 1. 

'^ Objects from lFbntana et al.l 1120101) are taken from their Table 1. For the objects of Vanzella ct al. (20ni), it is calculated from the Y- 



band magnitude corrected for Lya contribution and IGM attenuation, Vcont. For the objects of iSchenker et al.l 1 120111) and IPentericci et al 

(120111 '). it is calculated from the continuum magnitude rricont and redshift. The magniflcation factors are considered for the lSchenker et al 

pOlU ) objects. 

"^ For the objects of lFontana et aLl I I2010I) . the upper limit is estimated from t he Ict limiting flux at 94 85A, corresponding to z r^ 6.8 Lya 

around which the redshift distribution of their z-dropouts peaks (Figure 1 in ICastellano et ani2010al ) . 

+ The photometric redshift estimated by ISchenker et al.l 1 120111) or IMcLure et al.l 1 120111) . 

+ These are estimated bv lBradlev et aO H2011I V 

* The EW detection limits are well below EW= 25 A (IPentericci et al.|[20Til) . 
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l)iy 
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L } 
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Fig. 8. — Evolution in the fraction of strong LAEs in LBGs with -21.75 < A/uv < -20.25 (top panels) and -20.25 < Muv < -18.75 
(bottom panels) over 4 < 2 < 7. The left panels show the fraction of galaxies with EW larger than 25A, while the right panels show the 
fraction of those with EW larger than 55A. The filled square is our results, the open square is the results of IFontana et aLl I I2010I ) and 



^^^^__^^^^ ^_^^ _ ^^^^_^^^^^ ^_^ ^ pp e 

ISchenker et al.l ||20111 1. the cross is from 'VanzcUa ct al.' 1*2 01 ll 'l and IPentericci et al.l 1)20111 '). and the filled circle is the compos ite re sults. 
The filled diamonds are the results of Stark et al. (2011), and open triangle is the composite result of [Dow-Hveclund ct al.| l|2007l ^ and 
IStanwav et al.l 11200 71 V The filled square, open square, cross, and open triangle are shifted in redshift for clarity. The shaded area is derived 
by extrapolating the trend seen in lower redshifts to z r^ 7 l lStark et aLIPOilT) . 
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TABLE 6 
Summary of the samples 





EW^^" > 25A 


EWl;^" > 55A 


-21.75 < Muv < -20.25 


This Study 

Fontana et al. C2010') 

Vanzellaet al. (2011) 

Schenker et al. (2011) 

Pentericci et al. (2011) 


2/6T 
0/6 

2/2 
0/2 
0/7 


0/10 
0/6 
1/2 
0/2 
0/7 


-20.25 < A/uv < -18.75 


Schenker et al. (2011) 
Pentericci et al. (2011) 


0/1 
3/12 
0/3 


0/1 
2/12 
0/3 



Note. — For the sample of lSchenker et al.l II2011I ). ob- 
jects at 6.3 < 2 < 7.3 are considered. 

t In our sample, four objects have Lya EW limits larger 
than 25A. 



1.00 



I ' ' ' ' I 

_ v^i„<,=200km s-^ 
- v„i„d=25km s ' 








50 



100 



150 



EWo'^" 



Fig. 9. — Cumulative distribution function (CDF) of rcst-framc 
Lya EW for 2 ~ 7 galaxies. The filled circles are the composite 
re sults. The solid lines show the 2 = 7 CDF for the wind model 
of lDiikstra et al.l II2011I) with column density A^hi = 10^° cm'^^ 
wind velocity Vwind = 200 km s~^, and neutral hydrogen fraction 
XHi of 0.13, 0.41, 0.60, 0.80, 0.91 from top to bottom. The dashed 
lines are the same as the solid lines, except that fwind = 25 km 






(22 ±4%, 22 ±7%, 27± 
lations to z 

For an i ndepe ndent check 
IStark et al.l (I201TI) 



67 ±10% and X{^^^ 



and their cxtrapo- 
27±13%. 
o n the values from 
_ ..we refer t o IDow-Hygelund et al.l 
(|2007[) and Stanwa v et al.l (|2007D . who presented the re- 
sults of spe ctroscopy for i -dropout galaxies at z ■^ 6. 
IDow-Hygelund et ali ()2007D reported a total of 12 LBGs 
with bright UV continuum (-21.75 < Myv < -20.250- 
Three (one) of these galaxies are found to hayc Lya 
EWs larger than 25A (55 A). In an independent study, 
IStanway et al.l (j2007[) found one such bright galaxy (ID 
1042) with Lya EW of 23A. Combining these results, 
we estimate a Lya fraction of -^^lto = 23 ± 15% (3/13) 

20 BD 00, BD 03, BD 27, BD 44, BD 46, BD 58, BD 66, GOODS 
16 0, UDF PFs iO, UDF PFs 14, UDF PFs IDROPl, UDF PFsl 19 
(Table 4 of lDow-Hvgelund et al.1[2007l) . 



and X^^^ = 8 ± 8% (1/13) at z - 6, consistent with 
IStark et"aLl(|2011| ). Since they would not c h ange so much 
the extrapolations derived by IStark et all (|2011| ) due to 
the large statistical uncertainties of the combined results, 
we do not calculate extrapola tions using the com bined 
r esults as well as the re sults of IStark et al.l ()2011l) . 

iFontana et al.l ()2010D reported spectroscopic observa- 
tion of seven z-dropout galaxies, with six having bright 
UV continua (-21.75 < Muv < -20.25), and the other 
one having faint UV continuum (—20.25 < Muv < 
— 18.75). One of their galaxies shows Lya emission with 
EWo = 13A. Furthermore, IVanzella et~all (1201 ID per- 
formed spectroscopic observations of two z-dropout can- 
didates. Both of them have Muv < —20.25 and the 
two galaxies show L y a em ission with EWs of 64A and 
50A. iSchenker et al.l ()2011h showed the results of spec- 
troscopic observations of 19 galaxy candidates at 6.3 < 
z < 8.8 selected by t heir photometric rcdshift technique 
(jMcLure et al.|[20lTI ). They reported three Lya detec- 
tions with EWs of 65A, 65A, and 30A. iPentericci et al.l 
(|2011h presented the results of their spectroscopy for 11 
z-dropout candidates, among which one galaxy is consid- 
ered to be an interloper because of its multiple line de- 
tection. Seven galaxies have bright UV continua, while 
the remaining three galaxies are faint. Although two of 
the 10 candidates show Lya detection, their Lya EWs 
are lower than 25A. Note that Lya EW de tection lim- 
its for the objects of IPentericci et aTl (|2011[ ) is well be- 
low EWq^" = 25 A. Table [5] summarize s the spectro- 
scopic r esults of z-dro p out g al axies bv IFontana et al.l 
(l2010h. IVanzella et all (I20TI. ISchenker et al.l (I2011D 



and IPentericci et all (|2011l) . ISchenker et al.l (|201lD esti- 
mated the Lya fraction b ased on their results and the 
results of IFontana et al.l (|2010( ). correcting the effects 
of instrumentally-limited Lya visibilities in wavelength 
and strong OH airglow lines. They obtained X^^ = 

8.6j;?j% for UV-bright galaxies and X^^^ = 24 ± 14% 
for UV- faint galaxies. We estimate the Lya fra ction 
of z-dropout galaxies from IVanzella et al.l ()201lD and 
IPentericci et"al] (I201TI) to be X^^^ = 22ti?% (2/9) and 



ll_n% (1/9) for UV-bright galaxies at z - 7. 



v55 

For UV- faint galaxies, the upper limit of the Lya fraction 
is A:25^ < 33% (0/3) and X^^^ < 33% (0/3). Comparing 
estimates of the Lya fraction for z ~ 7 galaxies here, we 
find a spread due to small number statistics and possible 
field-to-field variation. Therefore, to minimize the Pois- 
son noise and cosmic variance, we combine estimates of 
Lya fraction from these independent studies. Combining 
these estimates assumes that the ficld-to-ficld variation 
is due to cosmic variance, and is not indicative of patchy 
reionization. 

We now combi ne our results fo r z ^ 7 galaxies here 
with th ose from IFontana et al. (|2010l ). IVanzella et aLl 

ISchcnkcr ct al] 



201 1|) . IPentericci et all (120111 ). and [Schcnkcr ct a! 



20111 ). to measure the Lya fraction using the spectra 



of z-dropout galax ies obtained by t hree d ifferent groups. 
For the objects of lSchenker et al.l (|2011| ). we use 14 ob- 
jects whose redshifts are in the range of 6.3 < z < 7.3, 
and take into account their correction factors for Lya vis- 
ibilities. Using the combined data, we calculate number- 
weighted mean Lya fractions: X'^^ = 17l9°% ^'^'^ 
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X^L. = 4 ± 4% for UV-bright galaxies, and X^^ 



Lya 



Lya 



We find that the Lya fractions X^^^ drop from 
to 7 based on the combined results. These fi nding s 
similar to those repor ted by iSchenker at all ()2011|) and 



24 ± 12% and Xl^^ = 16 ± 9% for UV-faint galaxies. In 
Figure [51 we present the evolution of Lya fraction with 
redshift for LBGs over the range 4 < z < 7. We also 
show, as shaded area, an extrapola tion of lower-z Ly a 
fraction to z ^ 7, which is derived bv lStark et al.l ()2011[ ). 

z - 6 
to 7 based on the combined results. These findings are 
similar to those repor ts 
iPentericci et all (120111) . 

IStark et al.l(|2010( ) reported that Lya fraction increases 
with redshift over 3 < z < 6 at fixed luminosity, and the 
trend is likely governed by redshift-dependent variations 
in dust obscuration, with additional contributions from 
kinematics and covering fraction of neutral hydrogen. In 
contrast to this increasing tendency of Lya fraction from 
z ~ 3 to 6, we have found a significant drop of Lya 
fraction from z ^ 6 to 7. This significant drop would 
suggest that Lya emission lines from z -^ 7 galaxies are 
scattered by the IGM whose neutral hydrogen fraction is 
higher than that at z ~ 6. We would witness the increase 
of neutral hydrogen fraction toward z ^ 7 in the cosmic 
reionization history. 

In Figure [U we find that X^^^ drops more strongly 
in the faint (Mi 



uv 



-19.5) galaxies than in the bright 
(Muv - -21) galaxies; X^sjz = 7;obs)/X§Jz = 
7; exp) = 0.36 ± 0.18, and 0.53 ± 0.33, respectively, 
where Xl^^{z = 7; obs) is the observed Lya frac- 
tion at z = 7, and X^^^(z = 7; exp) is the ex- 
pected Lya fraction at z = 7 derived by extrapolat- 
ing the trend seen in lower redshifts to z ^ 7. This 
magnitude dependence of -^Lya evolution could be ex- 
plained by different halo masses of galaxies and the 
surrounding IGM. Given that the clustering strength 
of drop out galaxies increases with their UV luminos- 
ity fe.g..lGiavalisco fc Dickinsonll2001l:lOuchi et al.l l2004l: 
lAdelberger et al.l 120051: iLee et al.ll2006[) . our results im- 
ply that the ionizing state of the IGM around galaxies 
hosted by less-massive dark matter halos changes later 
than that around galaxies hosted by massive dark mat- 
ter halos. This would suggest that reionization proceeds 
from high- to low-dens i ty environments (insi d e-out; e.g., 
Ciardi fc Ma"daul [20031: iSokasian et all [200l llliev et alj 



20061 c.f.. iFinlator et all 120090 rather t han from low- 



to hig h-density regions (ou tside-in; e.g.. iGnedinI 120001 : 
iMiralda-Escude et al. 2000l ). 

We compare our compo site results wit h those of model 
predictions derived by iDiikstra et al.l (|2011| ). which 
quantify the probability distribution function (PDF) of 
the fraction of Lya photons transmitted through the 
IGM, by combining galactic outflow models with large- 
scale seminumeric simulations of reionization. They as- 
sume that the IGM at z = 6 was fully transparent to 
Lya photons, and that the observed PDF for EWg^" at 
z = 7 is different only because of evolution of the ioniza- 
tion state of the IGM. Figure [HI compares their models 
with our composite results. Our results can be explained 
by an evolution of the neutral hydrogen fraction xhi be- 
tween z = 6 and 7; xhi is roug hly 0.6—0.9 at z ^ 7, w hich 
is similar to those repo rted bv ISchenker et al.l (|2011l ) and 
IPentericci et al.l (|2011| ). 

The above discussion assumes that the dropout sam- 



ples at different redshifts have similarly low contami- 
nation fractions from interlopers. We should, however, 
keep in mind that this assumption has not yet been jus- 
tified well, although for the z = 4 — 6 UV-faint sam- 
ples, IStark et al.l ( 20101 ) estimate d the contamination 
rate to be only 2 — 5 percent (c.f.. iStanwav et 311120081 : 
[Douglas et al.|[2009l ) . Another possible source of system- 
atic errors in our Lya fraction analysis is inhomogeneities 
among the dropout samples in the spectroscopic detec- 
tion limit and in the quality of the photometry used to 
estimate UV continua, both of which are difficult to fully 
take into account in our analysis. To reduce such inhomo- 
geneities, desirable is a systematic spectroscopic survey 
over an entire redshift range combined with deep pho- 
tometry redward of Lya wavelength. 

6. CONCLUSIONS 

In this paper, we have presented Keck/DEIMOS spec- 
troscopic observations of 11 z-dropout galaxies found in 
the SDF and GOODS-N fields. An emission line has been 
detected at 9500— lOOOOAin the spectra of three objects, 
one of which is the previously reported Lya emitter lOK- 
1 at z = 6.96. Since all the detected lines are singlet with 
a large positive weighted skewness, we have concluded 
that the three objects are Lya-emitting z-dropout galax- 
7.213, 6.965, and 6.844. Their Lya fluxes 



ies at z, 



spec 



and rest-frame Lya EWs are 2.5 x 10 

and 33A, 2 
-1 



-17 —1 —9 

' erg s cm 
X 10"" erg s-i cm-2 and 43A, 2.7 x lO'^'^ 
erg s~^ cm^^ and 43A, respectively. It should be noted 
that the z = 7.213 galaxy was confirmed by observations 
in two independent DEIMOS runs in 2010 and 2011 with 
three different spectroscopic configurations. This galaxy 
is detected in HST/WFCS F140W imaging as well as in 
Spitzer/IRAC 3.6 and 4.5 /im imaging. Stellar popula- 
tion modeling indicates that this galaxy has a very young 
age, a stellar mass < IO^Mq, and a high star formation 
rate of 30 — lOOM© yr~^, with strong nebular emission 
contributing to the fluxes in the IRAC bands. 

We then have measured the fraction of Lya-emitting 
galaxies at z ~ 7. To reduce statistical uncertain- 
ties and possible effects of field-to-field variance, we 
have combined our results with pr e vious z-dropout spec- 



trosc opic studies ([Fontana et al.l 20101: iVanzella et al.l 
2011!) including very recent ones ([Schenker et al.ll2011l : 
Pentericci et al.l [20 111 ) . We have obtained the z '^ 7 



Lya fraction of X^^ = 17 ± 10%, and X^^ = 4 ± 4% 



Lya 



for UV-bright galaxies (-21.75 < Muv < -20.25, or 
= 24±12%, and Xl^^ = 16±9% 



M, 



UV 



-21),andX§„ 



for UV-faint galaxies (-20.25 < Muv < -18.75, or 
Muv — —19.5). These low values indicate that the frac- 
tion of Lya-emitting galaxies drops from z ~ 6 to 7, in 
contrast to the reported increasing trend from z ~^ 4 to 
6. We have also found that Xl^^ drops more strongly 
in UV-faint galaxies than in UV-bright galaxies. These 
findings would suggest that the neutral fraction of the 
IGM significantly increases from z ^ 6 to 7, and that 
the increase is stronger around galaxies with fainter UV 
luminosities, which is consistent with inside-out reion- 
ization models where reionization proceeds from high- to 
low-density environments. 
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